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We study the properties of 2D fibre clusters and networks 
formed by deposition processes. We first examine the growth 
and scaling properties of single clusters. We then consider 
a network of such clusters, whose spatial distribution obeys 
some efi'ective pair distribution function. In particular, we 
derive an expression for the two-point density autocorrelation 
function of the network, which includes the internal structure 
of a cluster and the effective cluster-cluster pair distribution 
function. This formula can be applied to obtain information 
about nontrivial correlations in fibre networks. 
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There are many phenomena in nature that can be 
viewed as deposition problems, brought about through 
various transport mechanisms that bring particles to a 
surface. They include a multitude of processes such as 
deposition of colloidal, polymer and fibre particles 
Many deposition phenomena involve particles whose size 
is large compared to the their mutual interaction range, 
and so the main deposition mechanism is due to parti- 
cle exclusion. Among the most studied in this class is 
the Random Sequential Adsorption model ||l|,^. There 
particles are deposited on a surface and either stick or 
are rejected according to certain exclusion rules, with a 
maximum coverage (the "jamming limit") less than unity. 
This is in contrast to multilayer growth QJ^,^]. 

Processes of particle deposition may often begin from 
colloidal suspensions; solid particles suspended in a fluid. 
For some such systems, the interparticle repulsion is 
strong enough to prevent multilayer growth How- 
ever, the existence of dispersion forces can cause the par- 
ticles to flocculate, or aggregate, and to precipitate out 
of the suspension . For larger particles or clusters of 
particles, gravity often induces sedimentation out of the 
suspension ||]. 

A particularly interesting class of deposition-related 
problems which has received little attention is the depo- 
sition of fibres, fibre clusters, and the formation of fibre 
networks. Perhaps the most practical application of fibre 
deposition is that of paper-making. During its formation 
paper undergoes several stages, beginning as a colloidal 
suspension of fibres and ending up as a deposition of fi- 
bres on a surface, when the fluid is drained out. There 
have been many attempts to model the structure of fl- 



bre networks P,p|pO|-p^ . Much of this work has focused 
on the calculation of power spectra of ideal random fibre 
networks and their subsequent comparison with mass dis- 
tribution data obtained from paper-making experiments. 

In some sedimentation problems, such as in the mak- 
ing of laboratory paper sheets, the fibres can fiocculate 
while still in solution . This may occur for a variety of 
reasons, ranging from mechanical fibre-fibre interactions 
to hydrodynamic forces in the suspension ||] . In such cir- 
cumstances the resulting network may consist of clusters 
of fibres. Similar clusters are also expected to occur in 
other deposition problems. This is one of the main moti- 
vations for our work where we present a detailed study of 
the statistical properties of deposition processes involving 
individual fibres and clusters. 

We start with a simple model of deposition of fibres of 
length A and width lu on a 2D plane. Deposition begins 
from an initial seed fibre. At each deposition event only a 
fibre that overlaps at least one fibre already in the cluster, 
is kept. This growth rule is motivated by the adhesive 
sticking of fibres. The process forms a connected cluster 
of N fibres which is statistically spherically symmetric. 
Thus it is possible to define an average maximum radius 
R{N). Fig. |l| shows one configuration of a computer 
generated iV-cluster on a lattice, where N = 2000. We 
find that for >> 1 the radius satisfies 



R{N) = BNf^, 



(1) 



where the exponent (3 k, 1/3, over the range of N ex- 
amined and for all fibre geometries simulated, while the 
constant B depends only on fibre dimensions. Fig. ^ 
shows a plot of R{N) for two different types of fibre ge- 
ometries; A X w = 20 X 1 and A x w = 50 x 1 showing 
the corresponding plots on a log-log scale. We note that 
Eq. (|l|) also gives the number of fibres in a cluster directly 
from the radius. 

The exponent 1/3 in Eq. (|^) also arises from the 
following simple argument.]^ Assume we have a clus- 
ter of average radius R{N) composed of N fibres of lin- 
ear dimension a — V Aw. Next if we deposit an addi- 



^This result has been also independently derived by J. 
Astrom (J. Astrom, Pro Gradu avhandling, Abo Akademi 
(1989) (unpublished). 



1 



tional AN fibres onto the cluster its radius grows by 
AR and area by AA. The change in the cluster area 
AA = ia2AAf(27ri?a)/(7ri?2) ^ 2'kRAR, which in the 

continuum limit gives Eq. (|), withS = (3/27r)i/3 a. The 
factor 27ri?a/7ri?^ comes from the "active" zone at the 
edge of the cluster which is the only region where added 
fibres increase the area. 

Another measure to characterize the iV-cluster is the 
radial mass probability per unit length. This quantity is 
given by 



PN{r) = 



mN{r) 



(2) 



where mAr(r) represents the mass per unit area at a radial 
distance r from the center of the cluster. For large N, 
the normalized average mass probability density Pn (r) = 
(pAr(r)) can be written as 



PN{r) 



^R{N) 



f{r/R{N)), 



(3) 



where 7 = Jg°° f{u)du and where the average implies 
an ensemble average over all A'^-cluster configurations. 
Figure ^ shows a plot of Pn{'>') for five different values of 
N while the scaling function /(r) is shown in the inset. 
The scaling function / is universal for the values of A and 
uj studied here. The mass probability distributions have 
a peak-like component at the origin, continuing almost 
linearly for the greater part of their width. 

At the edge of the cluster, there is a fast decaying 
region that also scales as shown in Fig. ^ This is asso- 
ciated with the rough edge of the cluster which we can 
examine in the context of kinetic roughening of grow- 
ing interfaces . Defining fluctuations of the radius as 
W{N) = {{R{N)- R{N)yy/'^ , where R{N) is the radius 
for one cluster, we find it to grow asW^ N^, where x is 
consistent with 1/3 as expected since W/R = const. The 
roughness of the cluster is expected to alter the perco- 
lation threshold of a random fibre-cluster network from 
that of a distribution of solid disks [p^ . 

We now consider a disordered fibre network constructed 
by deposition of A^-clusters. This model is completely 
general and the only assumption is that there exists an 
effective spatial pair distribution function between any 
two clusters. One plausible application of such a network 
is the description of sedimentation of fibre suspensions 
containing floes of fibres, as discussed above. 

Starting from a distribution of ric fibre clusters, labeled 
by an index i, and where the i*'' cluster contains N{i) 
fibres, the mass distribution of the resulting network is 
M{x) = X]r=i ™JV(i)(-^ ~ where m^(j) is the mass 
density and Xi the coordinate of the center of the i*^ 
cluster. The two-point density autocorrelation function 
of this network is given by 



G(r) 



, ( (lAf(fc)l^ 



hm . 
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(4) 



where A is the area over which the 2D network is de- 
posited, is the inverse Fourier transform, and the 
average is taken over all network configurations. The 
Fourier transform of M{x) is given by 



N(m) ' 



(5) 



n—l m—1 



where /jv(n) is the Fourier transform of m7v(„) {x) evalu- 
ated over the area of the A^(n)-cluster. The coordinate 
in /jv(n) is defined relative to the center of the n"* clus- 
ter. Taking the ensemble average of Eq. (|^) we make the 
assumption that the description of individual clusters is 
independent of the positions of their centers. Thus 



-27rik-{xn—Xm) 



>x(/ 



(6) 



Next, making the assumption that internal structure of 
any given cluster is independent of the others, we write 
{lN{n)lNim)) = N {n)) (m)) noting that (/jv(„)) im- 
plies an ensemble average over all TOjv(n)(2^). This aver- 
age can be split into two parts: over all configurations 
where the n*'' cluster contains N{n) fibres and over the 
number N{n) itself. Breaking up the average in this way 
we define fhf^{r) = {m]y(^n){x)) to he the average mass 
density of an TV = {N{n)) cluster. Using Eq. (^) and not- 
ing that the total average mass in an A^-cluster is Nw / , 
where Wf is defined as the mass per fibre, we can show 
that 
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(r) = ^R{N)f{r/R{N)), 



(7) 



where 7 = 2 uf{u)du. Thus we have 

{lN(n)){I*N(rn)) = l^^ivd^DP: whcrc J^vd^l) is the Fourier 
transform of rhf^{r). To evaluate the exponential term on 
the right hand side of Eq. (^) we consider two cases: n = 
m and n ^ m. The first gives unity, while the second de- 
pends on the type of effective cluster-cluster interaction 
between clusters n and m. Defining Axn,m = Xn—Xm, we 
note that the ensemble average over all cluster-network 
realizations is equivalent to integrating e 
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over 



all Axn,m, weighted by a cluster-cluster pair probabil- 
ity density g(Axn,m)- We also assume, by isotropy, that 
g{Axn,m) is independent of the indices n and m, and ra- 
dially symmetric ijl^ . The average of the exponential on 
the right hand side of Eq. (|^) can thus be written as 



^ — 27Tik- Ax 
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(8) 



Combining the equations above we finally obtain 

|J^(fc)pg(fc), (9) 



(|M(fc)|2)=„,|J^(fc)|2 + 

ndric - 1) I 
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where g{k) is the Fourier transform of gd^D- Eq. (||) 
comprises two terms. The first represents the spectral 
density of Uc individual fibre clusters. Substituting this 
term into Eq. (Q) gives the two-point density correla- 
tion function of ric individual TV-clusters. The second 
term arises due to effective cluster-cluster interactions. 
The number of clusters Uc can be expressed in terms of 
a parameter 77, defined as the number of fibres per unit 
area p^ . In particular, we have ric — rjA/N. Writing 
N in terms of the average cluster radius and substitut- 
ing Eq. (||) into Eq. the following G'^(r) is finally 
obtained: 

4 . r°° 
7 -B-^'? Jo 

+ — / U^(k)g(k)kJoi2'n:kr)dk, (10) 
7^ Jo 

where Tl{k) (a dimensionless quantity) is given by 

/>oo 

U{k) = / f{z)zJo{2nkR{N)z)dz (11) 
"'0 

and where Jq is the Bessel function of order zero. In 
obtaining Eq. ( pO[ ) we normalized G{r) by dividing by 
the square of -qwf [|lO|,^ and noted that both n(A:) and 
g{k) are radially symmetric. 

Equation (nm is our main result for the fibre network. 
The first component in it describes the density-density 
correlations of a uniformly random distribution (g(r) — 
0) of fibre clusters. For this case it is easy to show that 
Eq. (|l^) becomes G^(r) = R{N)T{r/R{N)) where r(r) 
is Gaussian for small r and rapidly goes to zero as r ^ 
1. This dependence on R{N) in Eq. ([l^), allows us 
to extract the typical fibre-cluster size directly from the 
correlation function. 

The second term in Eq. ( p^ ) arises due to an effective 
pair distribution function acting between clusters. This 
function can lead to nontrivial long-range correlations in 
the fibre network. Such correlations could arise from hy- 
drodynamics in the sedimentation of particles from sus- 
pensions 1^ . We have examined a model pair distribution 
given by g{k) oc k^^. Figure ^ shows the second term of 
G7v(r) for x = 1-1) 1-3, 1.6 and R{N) = 3, where distance 
is measured in units of B. For r > R{N) (where the first 
term of Gf^{r) becomes negligible), we find G^(r) ~ r~" , 
where the correlation exponent a < 1 decreases contin- 
uously as X is increased from to 2. Since the function 
n(fc) « for fc > 1/R{N), we expect any effective pair 
distribution satisfying g{k) ^ k~^ for small k, to give 
essentially the same power law correlations at large r. 

In conclusion, we have presented a model for the statis- 
tical properties of disordered 2D fibre networks formed by 
deposition processes. The radius of fibre clusters scales 
with the number of fibres. Similarly, the density pro- 
files of a cluster obey scaling laws, with a universal scal- 
ing function for the cases studied here. Starting from 
the scaling form, we have derived the two-point density- 
density correlation function for a fibre network due to 



deposition of clusters, whose spatial distribution is given. 
In the case of uniform, random deposition of fibre clusters 
correlations exist over the range of the mean cluster size. 
However, an effective cluster-cluster interaction can give 
rise to long-range power law correlations. In this way 
the final structure of deposited fiber networks may give 
us information about the effective interactions between 
fibers in a suspension. 

We expect that our G^(r) can be compared with ex- 
periments on sedimentation of fiber-like particles in sus- 
pensions. A prime example of such a process is the forma- 
tion of laboratory paper sheets, whose final mass density 
can be measured by the beta-radiogram technique [ pT| . 
Our calculations will be compared with long-range corre- 
lations observed in such radiogram data in an upcoming 
publication. 
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FIGURES 

FIG. 1. A fibre cluster made of 2000 fibres of dimension 

X X uj = 50 X 1. The cluster is constructed by depositing 
fibres randomly, keeping only those that form one connected 
fibre cluster, beginning from an initial fibre. 



FIG. 2. A plot of R{N) on a log-log scale. In the bottom 
curve A X a; = 20 X 1, while for the top curve A x w = 50 x 1. 
The solid lines are drawn for comparison and have a slope of 
1/3. 



FIG. 3. The radial mass probability per unit area, 
for N = 1500,2000,2500,3000 and 4000 fibres, with 
A X o) = 50 X 1. The inset shows the scaling function f(x) vs. 
X = r/R. 



FIG. 4. The correlation function of Eq. (10) with 
g{k) = k^^, where x ~ 1.1,1.3,1.6 from bottom to 
top. In the large r limit power law behaviour r~™, with 
a = 0.92, 0.71, 0.4, emerges. 
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